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ABSTRACT 


In  this  report,  we  describe  an  algorithm  for  correlating 
measurements  from  several  sensors.  This  is  a  problem  area  in 
multiple  sensor  tracking  in  a  dense  target  environment.  It  is 
shown  that  the  correlation  problem  is  similar  to  the  assignment 
problem  in  operation  research  with  assignment  penalties  being 
equal  to  the  sufficient  statistic  of  the  generalized  likelihood 
ratio  test.  An  example  is  given  to  illustrate  the  correlation 
performance  for  two  optical  sensors  tracking  ballistic 
trajectories . 
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1.  INTRODUCTION 


In  many  target  tracking  applications,  one  is  confronted  with 
the  problem  of  multiple  sensors  simultaneously  tracking  a  number 
of  targets.  This  situation  brings  up  at  least  three  problem 
areas:  (1)  the  identification  of  a  time  sequence  of  measurements 
associated  with  the  same  target,  (2)  the  identification  of 
measurements  from  several  sensors  associated  with  the  same  target, 
and  (3)  efficient  algorithms  for  processing  multiple  sensor  data 
after  both  (1)  and  (2)  have  been  completed. 

The  first  problem  above  is  often  referred  to  as  the  problem 
of  multiobject  tracking  and  has  received  considerable  attention 
in  recent  years*.  The  second  problem  can  be  viewed  as  a  special 
case  of  the  first;  it  does  have, however ,  its  own  characteristics 
and  can  be  treated  separately.  This  is  the  subject  of  this  report. 
The  problem  (3)  above  was  the  subject  of  Ref.  [3]  where 
Kalman  filter  configurations  for  multiple  sensor  data  processing 
were  studied  in  detail. 

In  this  report,  we  present  an  algorithm  for  identifying  the 
same  object  from  measurements  made  by  several  sensors.  It  can 
be  shown  that  the  sensor  to  sensor  measurement  correlation 
problem  is  the  same  as  the  classical  assignment  problem  in 
operation  research,  [4]  —  [5] -  The  penalties  or  perfr rmance 

*See  Refs.  [1]  and  [2]  for  a  survey  of  open  literatures. 
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ratings  used  for  association  can  be  shown  to  be  the  sufficient 
statistic  of  the  multiple  hypothesis  generalized  likelihood  ratio 
test.  For  Gaussian  random  vectors,  the  sufficient  statistic 
follows  the  chi-square  distribution.  This  property  can  be  used 
to  eliminate  those  measurements  which  are  highly  unlikely  to  be 
correlated  with  the  same  target.  All  these  derivations  will 
be  shown  in  the  next  section.  In  the  third  section,  we  give  an 
example  showing  the  performance  of  correlating  measurements  of 
two  optical  sensors  tracking  multiple  ballistic  targets. 
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2 .  A  CORRELATION  ALGORITHM 

2.1.  The  Generalized  Likelihood  Ratio  Test 

Let  denote  the  j-th  measurement  of  the  i-th  sensor, 

assuming  that  there  are  M  sensors  taking  measurements  on  the  same 
N  targets.  The  problem  where  the  target  sets  observed  by  each 
sensor  are  not  identical  will  be  discussed  later. 

We  note  that  the  term  "measurement"  should  be  inter¬ 
preted  in  a  more  general  sense.  For  example,  y^  may  represent 
a  state  estimate  of  the  j-th  target  from  the  i-th  sensor.  In 
another  example,  a  subset  of  y^  may  be  actual  measurements  and 
the  remaining  subset  are  predicted  measurements  from  target  state 
estimators.  This  second  case  is  the  same  as  the  track  main¬ 
tenance  problem  of  multiple  sensor  tracking  in  a  dense  target 
environment. 

We  assume  that  all  y^'s  are  in  the  same  coordinate 
system.  For  a  specific  application,  the  required  transformation 
for  satisfying  the  above  assumption  can  be  defined. 

Let  denote  the  hypothesis  that  the  same  target  has 

generated  measurements  {y..  ;  i=l,...,M).  There  are  therefore  a 

1Ji 

M 

total  of  N  possible  hypotheses.  Let  x  denote  the  true  measure¬ 
ment  vector  corresponding  to  then  the  generalized  likelihood 
ratio  test  requires  the  following  decision  function  [ 6 ] 

d(y..  ;  i=l , . . . ,M) 

J  i 


(2.1) 


where  the  above  decision  function  is  obtained  by  using  the  maxi¬ 
mum  likelihood  estimate  of  x,  i.e.,  the  k  satisfying 


max  (p (y..  ;  i=l , . . . ,M/x) ) 
x 


(2.2) 


Substituting  k  in  the  density  p(^.  •  ;  i=l,...,M/x)  gives  the 

!  -j  i 

decision  function  (2.1).  Equations  (2.1)  and  (2.2)  form  the 
sufficient  statistic  for  correlating  measurements  from  multiple 
sensors. 


If  the  measurement  noise  vector  is  independent  between 
sensors  and  follows  a  Gaussian  density  function  with  zero  mean 
and  covariance  j ,  then  equation  (2.1)  becomes 


where 


M 

i=l,... ,M)  =  c  exp  (-1/2  ^  ^  (y_ii  -g)TPi;.  _1  (y^  ^  -x)  ) 

i=L  1  1  i 

(2.3) 


(2.4) 


Notice  that  the  k  above  is  the  maximum  likelihood  estimate  of  x 

assuming  that  all  ’s  are  correlated.  When  this  is  true,  it 

•^i 

is  the  compressed  measurement  discussed  in  [3],  Equations  (2.3)- 
(2.4)  give  the  sufficient  statistic  for  the  Gaussian  measurement 
noise  case. 

In  order  to  further  illustrate  the  above  results,  let 
us  consider  the  problem  of  correlating  measurements  from  two 
sensors.  Let  i=l,...,N)  and  (z y  j=l,...,N)  denote  the 
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measurements  of  the  first  and  second  sensors,  respectively. 
Furthermore,  we  assume  that  the  measurement  noise  vectors  of 


y.  and  z.  are  Gaussian  with  zero  mean  and  covariance 
*-1  -J 

respectively.  Using  these  assumptions,  then  Eq.  (2. 


=  c  exp  {-1/2  [  (^i-ft)TPi~1  (yi~x)  +  (Zj-x)TZ  j 


where 


One  can  now  simplify  Eq.  (2.5)  by  substituting  (2.6) 
This  is  done  with  the  following  derivations: 

(1)  Zi-a 

=  {I-[P  _1+Z  _1]  P4"1}xi-Ipr1  +  En_1]  E-Tlzn 

J.  J  X  X  X  J  J  — J 

=  [P._1+E  .“1]  (E  ._1y.-E  .-1z.) 

i  i  i  i 

_!  -1 

-  tEjpi  +IJ  (Zi-zj) 

(2)  (zi-a)Tpi"1(zi-^) 

=  (y.-z.)T[E  .+P.r1P.  (E  .+P.  ]-1  (y.-z  .  ) 

*-1-3  1  i  I  D  i  — D 

(3)  Similarly,  one  obtains 

z.-5c  =  [P.E.  ^+1}  ‘'"(y.-z.) 

— D  —  l  j  *-i  — d 

(zj-^)TEj*1(z.-a) 

=  (y.-z . )T[P.+E .]-1E . [P.+E .  ]-1 (y.-z.) 

'l-i  -j  i  jJ  D  1  D  — D 


P .  and  E  .  , 
i  D 

)  becomes 

(z ,-x) ]  } 

(2.5) 

(2.6) 

in  (2.5). 


(2.7) 


(2.8) 
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Substituting  (2.7)  and  (2.8)  into  (2.5)  yields 


d(Zi'Zj) 

=  C  exp  {-1/2 (Xi“£j )T [p^+z j ] _1 (Yi~lj  ) )  (2.9) 

Notice  that  all  information  for  decision  is  contained  in  the 
exponent  of  d(y^,Zj).  It  is  therefore  the  sufficient  statistic 
for  the  generalized  likelihood  ratio  test,  i.e.. 


=  (y  . -z  .  )  T  [P  .  +£  .  ]  '''(y.-z.) 
13  ■*-!  -3  1  3  ■*-!  —3 


(2.10] 


Notice  that  l ^ .  can  be  interpreted  as  a  weighted  distance  measure 

between  y.  and  z..  For  the  case  that  y.  and  z.  are  indeed  from 
£1  -3  -3 

the  same  target,  £ ^ ^  follows  a  chi-square  distribution.  For  the 
case  that  the  above  situation  is  not  true,  ^ -z..  is  not  zero  mean, 
8,^  therefore  has  high  probability  of  attaining  values  larger 
than  that  of  a  chi-square  random  variable.  This  fact  forms  the 
basis  of  using  chi-square  statistics  for  choosing  a  threshold 
for  pre-screening 

Consider  the  case  that  one  would  like  to  decide  which 
^  is  correlated  with  a  particular  z^ .  For  a  fixed  j,  the 
generalized  multiple  hypothesis  testing  procedure  will  select 
the  satisfying 


min  „  - 

i=l . N  ‘ij  ;  for  a  J. 


12.11) 
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If  one  simply  repeats  the  above  procedure  for  all  j ,  then 
ambiguous  situations  may  arise.  This  is  because  the  same 
measurement  of  a  given  sensor  may  be  selected  as  correlated  with 
several  measurements  of  the  other  sensor  while  some  measurements 
may  be  declared  uncorrelated  completely.  This  situation  may  be 
possible  when  sensors  do  not  have  perfect  resolution,  e.g.,  a 
sensor  may  report  only  one  measurement  for  several  closely 
spaced  targets  while  the  second  sensor  has  already  resolved 
these  targets.  Situations  like  this  can  be  treated  by  consider¬ 
ing  the  property  of  the  Chi-square  distribution  and  this  subject 
will  be  discussed  in  Section  2.3.  For  the  case  that  unambiguous 
matches  between  y^  and  must  be  made  (i.e.,  the  completely 
resolved  case) ,  this  gives  rise  to  the  classical  assignment 
problem. 


2.2.  The  Assignment  Problem 


Consider  Table  2.1.  Entries  of  the  matrix  can  be  £. . 's 

il 

of  Eq .  (2.10).  In  the  classical  assignment  problem  of  operations 

research,  the  may  represent  the  penalty  (or  payoff)  of 

assigning  the  i-th  person  to  work  on  the  j-th  job.  The  optimum 
assignment  is  selected  as  these  entries  giving 

min  y\. . 

i/i  ii 

i/  j 


(2.12) 


with  the  constraint  that  each  column/row  can  be  assigned  to  a 
row/column  only  once. 
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TABLE  2.1 


ENTRIES  OF  AN  ASSIGNMENT  PROBLEM 


Measurements  of  1st  Sensor 
ill  *2  *3  .  *N 


=  the  sufficient  statistic  of  assuming  and 
being  correlated. 
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Notice  that  the  above  matrix  is  similar  to  the  threat 

distribution  matrix  discussed  in  [8]  except  that  the  threat/targo t 

correlation  does  not  have  the  unique  match  constraint. 

A  trivial  approach  to  the  above  optimization  problem 

is  to  enumerate  all  possible  combinations  and  choose  the  one  with 

minimum  sum.  This  results  in  N!  trials.  A  classical  method 

called  the  Hungarian  algorithm  [4]  which  requires  substantially 

fewer  operations  has  been  known  for  many  years  in  the  field  of 

combinatorial  programming.  A  particular  method  of  implementing 

the  Hungarian  algorithm  given  by  Munkre  [5]*  requires  at  most 
3  2 

(11N  +12N  +3lN)/6  operations  representing  a  substantial  saving 
for  large  N's.  The  details  of  the  Hungarian/Munkre  algorithm 
will  not  be  repeated  here. 

There  is  a  "suboptimal"  approach  to  the  assignment 
problem  called  the  row  and  column  elimination  method.  This 
algorithm  proceeds  as  follows. 

(1)  select  an  arbitrary  row  (or  column) ,  find  its 
minimum  entry, 

(2)  search  the  corresponding  column  (or  row)  con¬ 
taining  this  entry  for  a  minimum. 

(3)  repeat  (1)  and  (2)  until  this  entry  is  the 
minimum  for  both  the  row  and  the  column  con¬ 
taining  it.  This  row  and  column  are  declared 
to  be  correlated. 


*The  authors  are  indebted  to  Dr.  R.  B.  Holmes  for  pointing  out 
this  reference. 


(4)  eliminate  this  row  and  column  from  further 
consideration,  then  repeat  ( 1 )  —  ( 3 )  until  all 
rows  and  columns  have  been  eliminated. 

The  above  algorithm  is  "suboptimal"  because  it  does  not  always 
give  optimal  solution  for  (2.12).  In  the  sensor-to-sensor  cor¬ 
relation  problem,  if  the  target  density  is  sufficiently  "low" 
such  that  a  which  does  not  correspond  to  the  same  target  will 

attain  a  much  larger  value,  then  the  row  and  column  elimination 
method  frequently  gives  the  optimum  solution. 

The  above  discussions  illustrate  the  problem  of  cor¬ 
relating  measurements  from  two  sensors.  For  more  than  two 
sensors,  the  computational  load  increases  rapidly  although 
the  problem  is  conceptually  straightforward.  In  a  multiple 
sensor  problem,  one  first  computes  the  sufficient  statistic 
using  (2.1)  —  (2.2) .  Let  it  be  denoted  as  l . 

i  1 ,  i 

where  M  is  the  total  number  of  sensors.  One  can  now  visualize 
the  l .  .as  entries  of  a  "super  matrix"  with  dimension 

(N  x  N  x  N  x. . .x  N)  where  there  are  M  N's.  The  correlation 
solution  is  given  by  the  combination  of  mutually  independent 
entries  achieving  the  minimum  sum.  It  is  interesting  to  note 
that  a  corresponding  Hungarian/Munkre  algorithm  for  the  afore¬ 
mentioned  problem  does  not  seem  to  be  available  yet.  An  ex- 

M-l 

haustive  search  method  will  require  enumerating  (N!) 
possiblities  and  this  number  looks  prohibitive  even  for  a 
modest  range  of  N  and  M.  The  computational  requirement  for  the 


row  and  column  elimination  method  is  still  reasonable  but  this 
method  works  well  only  for  a  modest  target  density. 

Finally,  we  comment  that  the  solution  criterion  given 
by  (2.12)  represents  one  of  many  possibilities.  A  more  general 
performance  index  can  be  chosen  as 

V  |  £ .  .|p  ,  for  p  >  1  (2.13) 

““ *  i  i  J 

i,  j 

max  H .  .  (2.14) 

•  ■  i/3 

i/3 

Notice  that  the  above  criterion  corresponds  to  the  minimum  norm 
problem  in  Banach  space.  If  one  chooses  (2.14)  as  the  criterion 
for  minimization,  this  corresponds  to  the  minimum  maximum  entry 
problem.  The  criterion  (2.12)  chosen  for  this  report  cor¬ 
responds  to  a  minimum  error  probability. 

2.3.  The  Chi-Square  Thresholding  Technique 

The  discussion  in  the  above  two  sections  constitutes  an 
essential  portion  of  a  sensor-to-sensor  correlation  algorithm. 

It  does  not  however,  consider  situations  when  (1)  sensor  field  of 
views  (FOV)  do  not  completely  overlap  and  (2)  the  sensor  resolu¬ 
tion  capability  is  limited.  For  these  situations,  we  consider 
an  approach  based  upon  the  Chi-square  statistics. 
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As  we  noted  earlier  for  the  Gaussian  measurement  noise 
case,  if  ^  is  formed  by  the  same  target,  it  is  a  Chi-square 
random  variable.  A  Chi-square  random  variable  x  with  n-th 
degree-of -freedom  follows  the  following  density  function  17]: 


p(x)  = 


-  , .  n/  n 
2(2:')  o  v  x 


x(n-l)/2e-x/2.r 


12 . 15) 


where  A  is  a  normalization  factor.  With  (2.15)  one  can  compute 

the  cumulative  probability  of  x. 

When  is  formed  by  two  different  targets,  then 

y.-z.  has  a  non-zero  mean.  For  this  case  the  ?  .  .  is  proportional 
-j  ij  p  K 

to  the  square  of  the  mean  vector,  it  will  therefore  attain  a 
value  larger  than  that  of  a  Chi-square  variable.  When  one  forms 
the  £ ^ j  matrix  (Table  2.1),  one  can  first  disqualify  those  large 
£^j's  because  they  are  statistically  unlikely  to  be  correlated. 
This  procedure  will  result  in  a  partially  filled  matrix.  Some 
parts  of  the  matrix  may  become  disjoint  from  the  others.  Their 
minimum  solutions  are  independent  of  the  rest,  they  can  therefore 
be  treated  independently.  Resolving  several  smaller  matrices  is 
computationally  more  efficient  than  resolving  one  larger  matrix. 
Some  of  the  submatrices  may  not  be  square  and  may  not  be 
completely  filled.  The  Hungarian/Munkre  method  can  be 
trivially  modified  to  handle  this  case. 
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A  tarqet  which  appears  only  in  the  FOV  of  one  sensor 


will  result  in  larqe  values,  and  will  therefore  likely  be 

discarded  in  the  above  thresholdinq  process. 

When  tarqet  density  is  high  so  that  some  sensors  may 
not  be  able  to  resolve  some  closely  space  targets,  the  -  values 
will  be  small  even  for  mismatched  cases.  In  order  to  reduce 
the  probability  of  leakage,  one  may  want  to  accept  one  measure¬ 
ment  of  one  sensor  to  be  correlated  with  several  measurements  of 
the  other  sensor.  For  examples,  let  £ ^  denote  the  i-th  measure¬ 
ment  from  the  first  sensor  and  j-th  measurement  from  the  second 
sensor  which  have  been  declared  correlated  using  the  method  of 
Section  2.2.  One  then  searches  for  those  entries  of  i-th  row  and 
]-th  column  which  are  within  a  certain  value  of  ^  .  All  these 
entries  will  be  accepted  as  possible  correlations.  This  method 
will  increase  the  correct  correlation  rate  for  the  case  with 
imperfect  sensor  resolution  and  high  target  density,  it  will 
nevertheless  also  increase  the  false  alarm  rate.  A  later  numeri¬ 
cal  example  will  illustrate  this  fact. 

The  above  discussion  establishes  a  method  for  reject¬ 
ing  highly  unlikely  correlations  and  for  accepting  multiple 
correlations.  This  method  may  be  ad  hoc  but  an  exact  approach 
seems  to  be  difficult  to  obtain. 

We  note  that  for  the  case  of  more  than  two  sensors, 
the  sufficient  statistic  for  Gaussian  measurement  noise  case  is 
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also  a  Chi-square  random  variable  and  the  discussions  above  also 
apply.  There  may  be  cases  such  that  a  target  appears  in  some 
sensors  and  does  not  appear  in  others.  Such  a  target  may  also 
obtain  a  small  sufficient  statistic  when  it  is  computed  with 
a  nearby  target.  This  is  because  the  sufficient  statistic 
(Eqs.  (2. 3) -(2. 4))  is  obtained  as  the  averaged  weighted  distance 
To  circumvent  this  problem,  one  may  first  want  to  perform  cor¬ 
relation  tests  for  all  two-sensor  pairs  and  isolate  those  target 
which  do  not  appear  in  all  sensors.  Those  targets  which  are 
identified  to  be  common  to  the  same  sensors  will  be  processed 
together. 


2.4.  Algorithm  Summary 

We  now  summarize  the  above  discussions  to  give  an 
algorithm  for  multiple  sensor  measurement  correlation. 

(1)  Choose  two  sensors,  compute  the  sufficient  statistic, 

between  their  measurements  using  (2.1)  —  (2.2)  ,  or 
(2.9)  for  the  Gaussian  measurement  noise  case. 

(2)  For  a  desired  leakage  probability,  a  threshold  can  be 
chosen  using  (2.15).  Apply  the  threshold  to  eliminate 
those  measurements  which  are  unlikely  to  be  correlated. 
Identify  those  targets  which  appear  to  only  one  sensor. 

(3)  Repeat  (1)  and  (2)  for  all  sensor  pairs,  identify  measure¬ 
ments  that  are  common  to  the  same  subset  of  sensors . 

(4)  Using  one  such  subset  of  sensors,  compute  the  sufficient 

statistic,  £.  .  ,  using  their  measurements  with 

1 1 ' 1 2  '  '  ’  *  '  ' 

Eqs.  (2.1)-(2.2),  or  (2.3)  —  (2.4)  for  the  Gaussian 
measurement  noise  case. 
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(5)  For  a  desired  leakage  probability,  a  threshold  for 

o. .  .  can  be  chosen  using  (2.15).  Apply  the 

11' ' ' • 

threshold  to  all  R  .  .  ’s.  For  those  cases  that 

11"  '  *  #1M 

are  above  the  threshold  the  corresponding  measurements 
are  declared  uncorrelated.  This  may  again  result  in 
some  measurements  not  being  correlated  at  all. 

(6)  Apply  the  Hungarian/Munkre  type  algorithm  to  process 
the  submatrices  resulted  from  thresholding.  This  gives 
the  correlation  solution. 

(7)  Let  R .  .  denote  an  entry  representing  a  set  of 

1 1 ' 1 2  '  '  '  *  '  1M 

correlated  measurements.  Fix  i2,...,iM,  search  those 

l .  .  for  all  i's  that  are  within  a  threshold 

1 ' x2 ' *  *  * ' 1m 

value  of  H.  .  .  .  Those  entries  which  are  within 

11'12'  '  *  "1M 

the  threshold  define  a  set  of  multiple  (or  redundant) 
correlations.  This  step  is  repeated  for  all  i.’s  and 
all  correlations  accepted  in  step  (6)  .  -1 

(8)  Repeat  steps  (4)- (7)  to  process  all  sensor  subsets 
identified  by  steps  (1)— (3) . 

The  above  steps  define  an  algorithm  for  correlating 
multiple  sensor  measurements.  The  correlation  criterion  is 
the  minimum  sum  of  weighted  distances  among  measurements  from 
each  sensor.  For  Gaussian  measurement  noise  vectors,  the  Chi- 
square  statistic  is  applied  to  eliminate  highly  unlikely  cor¬ 
relations.  A  method  (step  (7))  for  accepting  multiple  cor¬ 
relations  for  reducing  the  leakage  probability  is  described. 

We  note  that  a  multi-dimensional  Hungarian/Munkre  type 
of  algorithm  is  not  yet  available.  The  row  and  column  elimina¬ 
tion  method  may  work  well  for  modest  target  densities.  One  may 
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also  perform  correlation  test  for  all  pairs  of  sensors,  in  a 
sequential  fashion.  The  last  two  approaches  are  computationally 
much  simpler  although  they  are  not  optimum. 
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3. 


A  NUMERICAL  EXAMPLE 


In  this  section,  we  consider  an  example  which  illustrates 
the  target  handover  problem  between  two  sensors.  Consider  the 
case  of  two  optical  sensors  tracking  a  complex  of  ballistic 
targets.  The  first  sensor  tracks  a  subset  of  targets  for 
100  seconds  then  predicts  their  state  vectors  1000  seconds 
later.  The  second  sensor  tracks  another  subset  of  targets  which 
partially  overlaps  that  of  the  first  sensor  also  for  100  seconds. 
The  end  of  track  time  of  the  second  sensor  is  the  same  as  the  end  of 
the  prediction  time  of  the  first  sensor.  The  correlation  of  esti¬ 
mates  from  two  sensors  constitutes  a  typical  handover  problem. 

Initially,  the  number  of  targets  tracked  by  each  sensor  is 
equal  to  twenty  with  fifteen  of  them  being  common  to  both  sensors. 
The  initial  spacing  of  these  targets  is  strictly  uniform  and  very 
large  (with  a  5  to  1  target  spacing  to  tracking  uncertainty  ratio 
projected  to  the  angle  domain)  so  that  incorrect  correlations 
among  them  are  highly  unlikely.  These  targets  are  moving  with 
nearly  the  same  velocity.  The  target  density  is  increased  by 
introducing  more  targets  around  the  original  targets  where  the 
new  targets  are  randomly  distributed  with  a  uniform  density  over 
about  five  percent  of  the  original  target  spacing  and  moving  with 
a  very  small  re'ative  velocity.  We  will  call  an  original  target 
with  new  targets  surrounding  it  as  a  cluster.  Each  sensor  therefore  tracks 
twenty  clusters  where  the  initial  cluster  size  is  one  target. 


The  probability  of  correct  correlation  (P  )  and  the  pro¬ 
bability  of  false  correlation  (P^)  are  defined  as  follows 

p  _  number  of  targets  correctly  correlated 
c  ~  number  of  targets 

p  _  _ number  of  incorrect  correlations 

r  total  number  of  correlations  declared 

A  covariance  analysis  technique  (the  Cramer-Rao  bound) 
was  applied  to  obtain  the  track  accuracy  using  optical  sensors, 

[9] .  The  simulated  state  estimate  is  obtained  by  corrupting 
true  states  with  random  numbers  satisfying  the  Cramer-Rao 
covariance  bound.  A  correlation  performance  using  Pc  and  Pf 
is  obtained  with  cne  hundred  repetitions  on  each  target. 

The  and  Pf  as  a  function  of  number  of  targets  in  a 
cluster  is  shown  in  Fig.  3.1.  we  note  that  a  Chi-square  thres¬ 
hold  of  30  was  used  for  these  results.  Redundant  matches  were 
not  allowed.  In  this  case  the  number  of  targets  is  equal  to  the 
number  of  declared  correlations,  so  Pc  and  Pf  sum  up  to  one. 

Notice  that  Pc  drops  off  rapidly  as  the  number  of  targets 
in  a  cluster  increases.  This  is  due  to  the  small  spacing  of  targets 
within  a  cluster.  If  one  allows  for  redundant  correlations, 
the  Pc  can  be  improved  substantially  but  there  will  be  an  increase 
of  Pf.  Figure  3.2  shows  such  a  trade-off.  We  note  that  the  re¬ 
dundant  correlation  is  possible  by  accepting  those  targets  whose 


is 
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Probability  of  Correct  Correlation  (Pc> 


Chi-squares  are  within  a  certain  range  of  that  of  the  optimum 
choice.  Results  of  Fig.  3.2  are  obtained  with  this  approach. 

The  dotted  line  denotes  the  case  of  not  allowing  redundant  cor¬ 
relations.  Notice  that  the  increase  in  for  maintaining  a  90% 

P  is  rather  small, 
c 
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Probability  of  Correct  Correlation  (PJ 


1. 


Number  of  Targets 
in  A  Cluster 


Fig.  3.2.  Operating  characteristics  for  sensor- to- sensor  correlations 
with  target  density  as  a  parameter. 
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4. 


SUMMARY 


In  this  report,  we  have  presented  an  algorithm  for  cor¬ 
relating  measurements  of  several  sensors  tracking  in  a  multiple 
target  environment.  It  is  shown  that  the  correlation  problem  is 
the  same  as  the  assignment  problem  of  operations  research 
with  the  assignment  penalties  being  equal  to  the  sufficient 
statistic  of  the  generalized  likelihood  ratio  test.  For  the 
case  of  Gaussian  measurement  noise,  a  Chi-square  thresholding 
technique  is  applied  to  deal  with  the  problem  of  high  target 
density,  limited  sensor  resolution,  and  incomplete  overlapping 
of  sensor  field  of  views.  An  example  is  given,  illustrating  the 
correlation  performance  and  the  trade-off  between  correct  cor¬ 
relation  and  false  correlation  when  redundant  correlations  are 
allowed. 
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